	/*Alternative Analysen zum Gleichberechtigungsitem*/
	use turkey, replace

	/*L�nderdummies*/
	desmat sortcountry, full


	lab def sortcountry 5 DEO,modify
	lab def sortcountry 6 DEW,modify

	forvalues l=1/29 {
		local hilfe: label sortcountry `l'
		ren _x_`l' `hilfe'
		}

	lab def sortcountry 5 "DE-O",modify
	lab def sortcountry 6 "DE-W",modify
	



	label var DEO "DE-O"
	label var DEW "DE-W"


	label var AT  "AT"
	label var BE	 "BE"
	label var BG	 "BG"
	label var CZ	 "CZ"
	label var DK	 "DK"
	label var EE	 "EE"
	label var ES	 "ES"
	label var FI	 "FI"
	label var FR	 "FR"
	label var GR	 "GR"
	label var HR	 "HR"
	label var HU	 "HU"
	label var IE	 "IE"
	label var IT	 "IT"
	label var LT	 "LT"
	label var LU	 "LU"
	label var LV	 "LV"
	label var MT	 "MT"
	label var NL	 "NL"
	label var PL	 "PL"
	label var PT	 "PT"
	label var RO	 "RO"
	label var SE	 "SE"
	label var SI	 "SI"
	label var SK	 "SK"
	label var TR	 "TR"
	label var UK	 "UK"


	/*Rechtsstaatlichkeit*/
	/*Gemeinsame Regression mit Huber-White Standardfehlern*/

	reg rechtsstaatlichkeit religion mitbild hochbild wohnortgroesse male dekade berufstaetig AT-UK ,cluster(sortcountry ) nocons
	
	/*29 separate Regressionen?*/
	set more off

	reg rechtsstaatlichkeit religion mitbild hochbild wohnortgroesse male dekade berufstaetig  if sortcountry ==1
	capture matrix drop koeff
	matrix koeff = e(b)
	
	forvalues s=2/29 {
		reg rechtsstaatlichkeit religion mitbild hochbild wohnortgroesse male dekade berufstaetig  if sortcountry ==`s'
		matrix koeff = koeff \ e(b)
		}

	matrix colnames  koeff = ereligion         emitbild        ehochbild  ewohnortgroesse            emale          edekade    eberufstaetig econs


	capture drop listland
	gen listland=_n in 1/29
	label val listland sortcountry


	/*Ordered Logit f�r Gleichberechtigung*/
	/*Erfordert gologit2 von ssc*/
	

	capture matrix drop koeff
	gologit2 gleichberechtigung religion mitbild hochbild wohnortgroesse male dekade berufstaetig if sortcountry==1 , npl
	matrix koeff = e(b)

	forvalues s=2/29 {
		gologit2 gleichberechtigung religion mitbild hochbild wohnortgroesse male dekade berufstaetig if sortcountry==`s' , npl
		matrix koeff = koeff \ e(b)
		}


	/*Wohnortgr��e rekodieren*/
	recode wohnortgroesse (1/4=1 "<20T") (5/6=2 "20T-100T") (7/8=3 ">100T"),gen(urban)
	desmat urban
	rename _x_1 stadt
	rename _x_2 grossstadt



	
	set scheme s1manual



	
	/*Grafik erforder slideplot von ssc */


	
	slideplot hbar gleichberechtigung ,percent pos(3) neg(1 ) by(sortcountry )  plotregion(style(none)) ytitle(Prozent)

		graph export gleichberechtigung.eps, replace fontface(Times)


	/*Alternative Sch�tzungen mit gologit2/mlogit*/

	/*Variablen etwas anschaulicher kodieren*/
	gen invreligion=8-religion
	gen female = ~male
	gen jugend=60-alter

	/*Beruf? Interaktion mit Geschlecht?*/

	gen femaleberuf=female*beruf

	char sortcountry[pzat] ind(25)

	/*Schweden als Referenz festlegen*/
	
	mlogit gleichberechtigung invreligion mitbild hochbild stadt grossstadt  jugend female beruf femaleberuf ,base(2) cluster(sortcountry)

	gen hilfe=e(sample)
	est store pooled
	est save pooled, replace	
	
	/*Vollst�ndiges Modell*/
	set matsize 800

	desmat @invreligion*sortcountry  mitbild*sortcountry  hochbild*sortcountry  stadt*sortcountry  grossstadt*sortcountry  @jugend*sortcountry female*beruf*sortcountry
	mlogit gleichberechtigung _x* ,base(2)


	
est store full
est save full, replace


	/*Modell mit l�nderspezifischen Konstanten*/
	
	mlogit gleichberechtigung invreligion mitbild hochbild stadt grossstadt  jugend female beruf femaleberuf AT-RO SI-UK,base(2)
	est store dummies
	est save dummies,replace
	

	/*Intercept-Only-Modell*/
	mlogit gleichberechtigung AT-RO SI-UK if hilfe,base(2)

	est store intercepts
	est save intercepts,replace
	

	/*Besser gologit2?*/
	gologit2 gleichberechtigung AT-RO SI-UK if hilfe

	est store ordintercepts
	est save ordintercepts,replace

		
	gologit2 gleichberechtigung invreligion mitbild hochbild stadt grossstadt  jugend female beruf femaleberuf AT-RO SI-UK
	est store dummies
	est save dummies,replace

	gologit2 gleichberechtigung invreligion mitbild hochbild stadt grossstadt  jugend female beruf femaleberuf , cluster(sortcountry)
	est store ordpooled
	est save ordpooled, replace

	/*Variablenlabel anbringen*/
	label var invreligion "-Religion"
	label var mitbild "Bildung: mittel"
	label var hochbild "Bildung: hoch"
	label var stadt "Wohnort $\geq$ \np{20000}"
	label var grossstadt "Wohnort $\geq$ \np{100000}"
	label var female "Frau"
	label var beruf "Beruf"
	label var jugend "60-Alter"

	label var gleichberechtigung "Jobs are scarce: giving men priority"

	
	mlogit gleichberechtigung invreligion mitbild hochbild stadt grossstadt female  jugend beruf AT-RO  SI-UK, cluster(sortcountry) base(2)

	fitstat

	local mfadj=r(r2_mfadj)
	local coxsnell=r(r2_ml)
	local nagelkerke=r(r2_cu)
	
		/*So ist es eigentlich h�bscher*/
	/*Vorsicht: die Pseudo-R2 sind zu hoch, weil hier der Vergleich mit dem Intercept-Only-Modell irref�hrend ist*/

	mlogit gleichberechtigung invreligion mitbild hochbild stadt grossstadt female  jugend beruf AT-RO SE SI-UK, cluster(sortcountry) base(2) nocons


	
	/*Achtung, spost scheint nur unter version 9.2 richtig zu funktionieren*/

	
version 9.2
	mlogit gleichberechtigung invreligion mitbild hochbild stadt grossstadt female  alter beruf AT-RO SE SI-UK, cluster(sortcountry) base(2) nocons
	
	/*Tuerkei, gebildete Gro�st�dterin, nicht fromm*/
	prgen alter,gen(turkfrau) from(18) to(80) x(female 1 beruf 1 hochbild 1 grossstadt 1 invreligion 1 TR 1) rest(min)
	prgen alter,gen(turkfrauhlg) from(18) to(80) x(female 1 beruf 1 hochbild 1 grossstadt 1 invreligion 7 TR 1) rest(min)

	/*Tuerkei, gebildeter Mann, nicht fromm*/

	prgen alter,gen(turkmann) from(18) to(80) x(beruf 1 hochbild 1 grossstadt 1 invreligion 1 TR 1) rest(min)
	prgen alter,gen(turkmanhlg) from(18) to(80) x(beruf 1 hochbild 1 grossstadt 1 invreligion 7 TR 1) rest(min)

	/*Religion vertauscht!*/
	lab var turkfraup3 "Frau, sehr religi�s"
	lab var turkfrauhlgp3 "Frau, nicht religi�s"
	lab var turkmannhlgp3 "Mann, nicht religi�s"
	lab var turkmannp3 "Mann, sehr religi�s"

	graph twoway line turkfraup3 turkfrauhlgp3 turkmannp3  turkmannhlgp3  turkmannhlgx,xtitle("Alter") ytitle("Pro Gleichberechtigung") plotregion(style(none)) lpattern(dash dash solid solid) lcolor(black gs13 black  gs13) note("T�rkei")
	graph rename turkey,replace
	
	graph export turkey1.eps, replace fontface(Times)

	/*Das gleiche noch mal fuer Schweden. In beiden F�llen fehlen ein paar drop-Kommandos f�r die x-Variablen*/

	prgen alter,gen(swedfrau) from(18) to(80) x(female 1 beruf 1 hochbild 1 grossstadt 1 invreligion 1 SE 1) rest(min)
	drop swedfx

	prgen alter,gen(swedfrauhlg) from(18) to(80) x(female 1 beruf 1 hochbild 1 grossstadt 1 invreligion 7 SE 1) rest(min)

	prgen alter,gen(swedmann) from(18) to(80) x(beruf 1 hochbild 1 grossstadt 1 invreligion 1 SE 1) rest(min)
	prgen alter,gen(swedmanhlg) from(18) to(80) x(beruf 1 hochbild 1 grossstadt 1 invreligion 7 SE 1) rest(min)


	lab var swedfp3 "Frau, sehr religi�s"
	lab var swedfrauhlgp3 "Frau, nicht religi�s"
	lab var swedmanhlgp3 "Mann, nicht religi�s"
	lab var swedmannp3 "Mann, sehr religi�s"

	graph twoway line swedfp3 swedfrauhlgp3 swedmannp3  swedmanhlgp3  swedmanhlgx,xtitle("Alter") ytitle("Pro Gleichberechtigung") plotregion(style(none)) lpattern(dash dash solid solid) lcolor(black gs13 black  gs13) note("Schweden")

	graph rename sweden, replace

	
	graph combine turkey sweden,cols(1) ycommon ysize(5) xsize(4)
	graph export turkey-sweden.eps,replace fontface(Times)
